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Abstract 

In this paper we study the steady state of the fluctuations of the surface for 
a model of surface growth with relaxation to any of its lower nearest neighbors 
(SRAM) [F. Family, J. Phys. A 19, L441 (1986)] in scale free networks. It is known 
that for Euclidean lattices this model belongs to the same universality class as 
the model of surface relaxation to the minimum (SRM). For the SRM model, it 
was found that for scale free networks with broadness A, the steady state of the 
fluctuations scales with the system size N as a constant for A > 3 and has a 
logarithmic divergence for A < 3 [Pastore y Piontti et al, Phys. Rev. E 76, 046117 
(2007)]. It was also shown [La Rocca et al, Phys. Rev. E 77, 046120 (2008)] that 
this logarithmic divergence is due to non-linear terms that arises from the topology 
of the network. In this paper we show that the fluctuations for the SRAM model 
scale as in the SRM model. We also derive analytically the evolution equation for 
this model for any kind of complex graphs and find that, as in the SRM model, 
non-linear terms appear due to the heterogeneity and the lack of symmetry of the 
network. In spite of that, the two models have the same scaling, but the SRM model 
is more efficient to synchronize systems. 

Key words: complex networks, interface growth models, transport in complex 
networks. 
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Recently, much effort has been devoted to the study of dynamics in complex 
networks. This is because many physical and dynamic processes use complex 
networks as substrates to propagate, such as epidemic spreading [I], traffic 
flow Pf5f4"] and synchronization |5|5] . In particular, synchronization problems 
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in networks are very important in many fields such as the brain network [7], 
networks of coupled populations in epidemic outbreaks [8] and the dynamics 
and fluctuations of task completion landscapes in causally-constrained queuing 
networks [SJ. Synchronization deals with the optimization of the fluctuations 
in the steady state of some scalar field h, that can represent the neuronal 
population activity in brain networks, infected population in epidemics and 
jobs or packets in queuing networks. It is particularity interesting to under- 
stand how to reduce the load excess in communication networks in the steady 
state. This problem can be mapped into a problem of non-equilibrium surface 
growth where h is a random scalar field on the nodes that could represent the 
total flow on the network and the load excess could represent the overload of 
flow that a node should handle. Then, a way to reduce the load excess is to 
reduce the fluctuations of that scalar field. Recently, Pastore y Piontti et. al 
[TU] used the model of surface relaxation to the minimum (SRM), that allow 
balancing the load, reducing the fluctuations in scale free (SF) networks with 
degree distribution given by P{k) ~ k~ x , (k > k min ) with k the degree of a 
node, k m in the minimum degree that a node can have, and A the broadness 
of the distribution [llj. Given a scalar field h on the nodes, that in surface 
problems represents the interface height at each node, the fluctuations are 
characterized by the average roughness W(t) of the interface at time t, given 

by W = W{t) = {l/NJ2?=i(hi- (/i» 2 } V \ where K = h^t) is the height 
of the node i at time t, (h) = (1/N) Y^Li hi, N is the system size, and {.} 
denotes average over configurations. 

The aim of this paper is to study the steady state of W = W s for the model of 
surface growth with relaxation to any of its lower nearest neighbors (SRAM) 
[T2] . We find that this model has the same behavior with the system size as the 
SRM model for every A, even though the SRM model is more efficient to reduce 
the fluctuations and to enhance synchronization than the SRAM model, as we 
show later. Moreover, we derive analytically the general evolution equation for 
the SRAM model for any kind of random graph. 



1 Surface relaxation models 



In the SRM model, at each time step a node i is chosen with probability 1/N. 
If we denote by v% the nearest-neighbor nodes of i, then 



if < 



(1) hi < hj, Vj G Vi 

(2) hj <h n Wn^ j e Vi 



hi = hi + 1 



hj = hj + 1. 
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The SRAM model has the same rule (1) as the SRM model, but the second 
rule is different: the chosen node can relax to any of its lower m neighbors 
with probability 1/m. Then, the rules for the SRAM model are 



if 



(1) hi < hj, Vj G Vi 



=>- hi = hi + 1, else 



(2) 3 m nodes j G Vi with hj < hi =>■ hj = h j + 1 with probability 1 / 



It is known that in Euclidean lattices of typical linear size L the SRM and the 
SRAM models belong to the same universality class [T21T3] . with 




where t s is the saturation time which scales as t s ~ L z . Here, the exponent (3 
is the growth exponent, a is the roughness exponent and z is the dynamical 
exponent that characterizes the growth correlations given by z = a/j3. For 1 
dimension, these exponents are (3 = 1/4, a = 1/2 and z = 2. Moreover, both 
growth models belong to the same universality class as the Edward- Wilkinson 
(EW) equation, 



where v is the coefficient of surface tension and r](x) is a Gaussian noise with 
zero mean and covariance given by 

{rj(x, t)rj(x, t')} = 2D5(x - x')5(t - t) . 

Here, D is the diffusion coefficient and is taken as a constant. The fact that 
these models are represented by the same phenomenological equation is due 
to the symmetry of those models on the underlying Euclidean substrate [T3] , 
The extension of the EW equation to any unweighted graph with iV nodes is 
described by 





dh(x, t) 
dt 



uV 2 h(x,t) +r}(x,t) , 



dhjjt) 
dt 



N 



= V 



^Ay(/i i (t)-/i i (t))+7 7i (t) 



(2) 



where i and j are nodes of the graph, {Aij} is the adjacency matrix (Aij = 1 
if i and j are connected and zero otherwise), v represents the same as in 



Euclidean lattices and f]i(t) is a white Gaussian noise with zero mean and 
covariance given by 



being D the same as in Euclidean lattices. 

In Ref. [10J it was found that the saturation regime of W s in SF networks 
scales with N as 



It was also shown that the EW equation given by Eq. ([2]) predicts that in the 
thermodynamic limit W s ~ l/(k) for any random graph [6J. Then, the un- 
weighted EW equation in random graphs cannot describe the SRM model. La 
Rocca et. al [15], using a temporal continuous approach, derived the evolution 
equation that describes the SRM model. They found that the logarithmic di- 
vergence for A < 3 cannot be explained by the unweighted EW equation [See 
Eq. in graphs. The equation derived in Ref [T5] contains non-linear terms 
and weights that appear as a consequence of the heterogeneous topology that 
a SF has for A < 3, even though the network is unweighted. The heterogene- 
ity breaks the symmetry h — > — h of Eq. (|2"1) . For A > 3 the heterogeneity is 
not strong enough and the non-linear terms are negligible for the system sizes 
studied there, and the behavior of the fluctuations becomes well described 
by a weighted EW equation. It is not unexpected that transport processes 
in random heterogeneous graphs behave differently than in Euclidean lattices 
due to the fact that the nodes with high degrees (hubs) play a mayor role in 
transport. For example, reaction-diffusion processes behave very differently in 
homogeneous lattices than in SF networks with 2 < A < 3 due to the pres- 
ence of hubs which control the behavior for long times [TBffTT] . The hubs are 
responsible of a super diffusive regime because they diminish the distances. 



2 Saturation results for the SRAM model 

We construct our networks using the Molloy-Reed (MR) algorithm [TB], with 
kmin — 3 in order to ensure that the network is fully connected. The initial 
conditions for the scalar field {h} were drawn from a random uniform distri- 
bution between [0, 1]. In Fig. [fl we plot W 2 as a function of t for A = 3.5 and 
2.5 and different values of N. In the insets of Fig. [1] we plot W s as a function 
of N. We can see that for A = 3.5, W s increases, but asymptotically goes to a 



{Vi{t)Vj(t')} = 2D SijSit-t) , 




const, for A > 3 



In iV for 2 < A < 3 . 



(3) 
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constant and all the iV dependence is due to finite-size effects as in the SRM 
model [10]. However, for A = 2.5 we find that W s ~ In AT, as in the SRM 
model. This is shown in the inset of Fig. UJ6) in log-linear scale. Then, W s 
for both models scales in the same way for SF networks, with a logarithmic 
divergence for A < 3 and as a constant for A > 3. 



3 Analytical Evolution Equation 



Next, we derive analytically the evolution equation for {h} of the SRAM model 
for any kind of random graphs. 



The procedure chosen here is the same as the one used in Ref. [15] and is 
based on a coarse-grained (CG) version of the discrete Langevin equations 
obtained from a Kramers-Moyal expansion of the master equation [T91ll3f2Uj . 
The discrete Langevin equation for the evolution of the height in any growth 
model is given by [T3f2T)] 



9hi 1 „ , . , 

— = -d + m , 4 
at r 

where Gi represents the deterministic growth rules that cause the evolution of 
the node i, r — NSt is the mean time to grow a layer of the interface, and r\i 
is a Gaussian noise with zero mean and covariance given by [13120] 

{Vi(t)vj(t')} = -G i d ij 6(t-t') . (5) 
r 



If k 1 represents the degree of node 7, we can write Gi more explicitly as 

N 

G i = uj i + Y d K 1 + + ... + J}') , (6) 
3=1 



where Ui is the growth contribution by deposition on the node i and uj 1 is 
the growth contribution to the node i by relaxation from its neighbor j with 
probability 1/m, being m the number of neighbors of the node j with smaller 
heights than the node j. Then, 

= n 1 1 - @ (^- - ^ > 

jevi 

Uj) = [1 - - hj)\ II [1 - ®( h 3 - h n)} , 



5 



u] = [1 - Q(h { - fy)] n [1 - ®( h i - M 5 [1 - e(/i m - hj)] 



[i - Q(hi - h^] n t 1 - e ( /i i - M 7T [i - ®(h m - hi)] [i - e(h t 



f = [i - efa - ^)] n [i - ec/in - ^-)] r 



Here, is the Heaviside function given by Q(x) — 1 if x > and zero 
otherwise, with x = h t — h s = Ah. Without loss of generality, we take r = 1 
and assume that the initial configuration of {hi} is random. A schematic plot 
of the growing rules are shown in Fig. [5] where the rule (1) represent Ui, (2) 
represent and (3) represent w|. 

In the CG version x — > 0; thus after expanding an analytical representation 
of Q(x) in Taylor series around x = to second order in x, we obtain 



N N 

d= a k * + J2 c a + ci^ 1 " 1 Y.Aijihj-hi) 

3=1 i=i 

Ci ^ ci N N 

+ — J2 c v( h j- h i) + — Y, T i3 J2 A jn (h n -hi) 



C-2 



N 



j=l n=l,n=£i 



— ^ Cij(hj — h, 



? - a^ 1 



i=i 



N 



X! Aij{hj — hi 



A/ 



+ E 



o oQij ^'ij 



2a 2 ^ a 13 



N 



n=l,n^i 
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a ki 2 c 2 
2 

c 2 N 



N 



7, ~ hi 



„2 AT 



2a 2 



A" 



7, Aj n (h n — foj) 



2 



7 ] Ajnihn hj 



(7) 



where a = (1 — c ), c , Ci and c 2 are the first three coefficients of the expansion 
of the Q(x) and 
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"J 



Q 



'j 



2(1-2*0 + tf + kj 



(kj - 2)(kj - l)kj 
are different "weights" on the link ij introduced by the dynamics. 

In our equation the non-linear terms in the difference of heights arise as a 
consequence of the lack of a geometrical direction and the heterogeneity of the 
underlying network. This result is very similar to the one found for the SRM 
model in SF unweighted networks [15], where the non- linear terms appear due 
to the heterogeneity of the network. 

For the noise correlation [See Eq. (j5J)], up to zero order in Ah [T3f2"0] we obtain 
{r h {t)r ]j {t')} = 2D{k i )8 ij 8{t-t') with 



D{k t ) = -[a* + 1 £C ij ]. (8) 

/ 3=1 

Notice that all the coefficients of the equation depend on the connectivity of 
node i, i.e., on the topology of the underlying network. This dependence on 
the topology is expressed as weights on the links of the unweighted underlying 
network that appears only due to the dynamics on the heterogeneous network. 



4 Numerical results of the analytical evolution equation 



We numerically integrate our evolution equation for SF networks taking into 
account the linear terms and the first non- linear term of Eq. (|7|). This is 
because when non-linear terms are considered, the numerical integration al- 
gorithms we use has numerical instabilities. This is still an open problem to 
be solved in the future. For all our integrations we used the Euler method 
with the representation of the Heaviside function given by Q(x) = {1 + 
tanh[[/(x + z)]}/2, where U is the width and z = 1/2 [20], and random ini- 
tial conditions. With our choice of the representation of the Heaviside func- 
tion, we obtain: c = [1 + tanh(£//2)]/2, d = [1 - tanh 2 ([//2)] U/2, and 
c 2 = [-tanh([//2) +tanh 3 (f//2)] U 2 /2. 

In Fig. [31 we plot W 2 as a function of t, obtained from the integration of 
Eq. flU with Eq. © and D(h) given by Eq. © for A = 3.5 and 2.5 and 
different values of N, with k min = 3. For the time step integration we chose 
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At <C i/k max according to Ref. [21], where k max ~ iV 1 /^ -1 ) is the degree 
cutoff for the MR construction. In the inset figures we plot W s as a function 
of N. We can see that for A = 3.5, W s increases, but asymptotically goes to 
a constant and all the N dependence is due to finite-size effects. However, for 
A = 2.5 we found a logarithmic divergence of W s with N, as shown in the 
inset of Fig. [3](6) in log-linear scale. The fit of W s with a logarithmic function 
for A = 2.5 shows the agreement between our results and those obtained for 
the SRAM model in SF networks for A < 3. Then, our equation reproduces 
correctly the behavior of W s for the model for any A > 2. Notice that for A > 3 
and the system sizes studied here, the non-linear terms do not contribute and 
the process can be described by a weighted EW equation. 



5 Discussions 



The behavior of the SRAM and SRM models in the steady state are the same 
and both evolution equations are in agreement with the fact that all the co- 
efficients depend on the connectivity of a node, i.e., on the topology of the 
underlying network, and the weights appear only due to the dynamics on the 
heterogeneous network. Another similarity between both models is that for 
A > 3 the non-linear terms do not play any role for the systems size studied 
here, then both process are well described by a weighted EW equation. How- 
ever, the equations for both models are different [TJ)] and the main difference 
appears in the weights that produce different values of W s among them. For 
the synchronization problem, the SRM model is more efficient than the SRAM 
model, as can be understood from a lower W s shown in Fig. H] where we plot 
W 2 as a function of t in log-log scale for the SRM and SRAM models for 
N = 1024 for (a) A = 3.5 and (6) A = 2.5. We can see that the SRM model 
reaches the saturation regime faster than the SRAM model and W s for the 
first one is much lower than W s of the latter. This means that for the same 
system size, the process that will be better for synchronizing is the SRM, since 
it has less fluctuations on its scalar fields. These observations can be explained 
as follow: in both models, the nodes with low degree control the process all 
the time because they are more abundant and make a major contribution to 
the growth of the hubs. Then, we expect that the growing contribution of the 
hubs will be by relaxation from their neighbors with lower degree because our 
networks are disassortative for A < 3 (due to the MR construction) [22]. This 
is because in the SRM model at the initial stages the hubs grow faster than 
in the SRAM model. As the hubs are more important in the SRM model than 
in the SRAM model, the height-height correlation length should growth faster 
allowing to reach the saturation regime earlier. Notice that in the SRM model 
the nodes relax always to the minimum, while in the SRAM model the relax- 
ation takes place at any randomly chosen neighbor with smaller height (not 
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necessarily the minimum) than the chosen node. Then, if we have to chose 
one of these models as a synchronization process, it is more efficient to use 
the SRM model. 



6 Conclusions 

In summary, we studied the SRAM model in SF unweighted networks and 
found that for A > 3, W s scales as a constant and the N dependence is due 
to finite-size effects, while for A < 3 there is a logarithmic divergence with 
N, the same as in the SRM model. Then, the SRAM and SRM models still 
scale in the same way for SF networks. We derived analytically the evolution 
equation for the SRAM model for any network and find that even when the 
underlying network is unweighted, the dynamic introduces weights on the 
links that depend on the topology of the network. This equation contains 
non-linear terms and considering the linear and only the first non-linear term 
in the integration of the evolution equation, we recover the scaling of W s with 
N for any A > 2. And last but not least, we found that even though the two 
models have the same scaling, for synchronization problems the SRM model 
is more efficient because it reaches the steady state faster than the SRAM 
model and its fluctuations are much lower. 
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Fig. 1. W 2 as a function of t for the SRAM model: (a) A = 3.5 for iV = 64 (O), 128 
(□), 256 (o), 512 (A), 1024 (y), 1536 (★) and 2048 (X). (b) A = 2.5 for N = 1024 
(O), 1280 (□), 1536 (o), 1792 (A), and 2048 (y)- In the inset figures we plot W s 
as a function of in symbols. The inset figure of (b) is in log-linear scale and the 
dashed line represents the logarithmic fitting of W s with N. 



1/2 1/2 
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Fig. 2. Schematic plot of the growing rules for the SRAM model in one-dimensional 
Euclidean lattice. The solid arrows indicate the contributions to the growth of site 
i due to deposition (1) and diffusion from the nearest neighbors (2) and (3). Notice 
that in the case (3) site i growths with probability 1/2. 
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t t 

Fig. 3. W 2 as a function of t for the integration of the evolution equation using the 
linear terms and the first non-linear term in Eq. ([7]): (a) A = 3.5 for N = 64 (0)> 
128 (□), 256 (o), 512 (A), 1024 (y) and 1536 (*). (b) A = 2.5 for N = 192 (O), 
256 (□), 384 (o), 512 (A), and 768 (y). In the inset figure we plot W s as a function 
of N in symbols. The inset figure of (6) is in log-linear scale and the dashed line 
represents the fitting with W s ~ luN. 
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Fig. 4. as a function of i in log-log scale for the SRM (O) and SRAM (□) 
models for N = 1024 for: (a) A = 3.5 and (b) A = 2.5. 



12 



